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The mixing mechanism of axial- vectors Dsi(2460) and -Dsi(2536) is studied via intermediate 
hadron loops, e.g. D* K, to which both states have strong couplings. By constructing the two-state 
mixing propagator matrix that respects the unitarity constraint and calculating the vertex coupling 
form factors in a chiral quark model, we can extract the masses, widths and mixing angles of the 
physical states. Two poles can be identified in the propagator matrix. One is at y^i = 2454.5 MeV 
corresponding to _Dsi(2460) and the other at y^s = (2544.9— l.Oi) MeV corresponding to Z)si(2536). 
For Dsi(2460), a large mixing angle 9 = 47.5° between "^Pi and ^Pi is obtained. It is driven by 
the real part of the mixing matrix element and corresponds to 9' — 12.3° between the j = 1/2 
and j = 3/2 state mixing in the heavy quark limit. For Z)ai(2536), a mixing angle 9 — 39.7° 
which corresponds to 0' = 4.4° in the heavy quark limit is found. An additional phase angle 
(j) = —6.9° ~ 6.9° is needed at the pole mass of Dai (2536) since the mixing matrix elements are 
complex numbers. Both the real and imaginary part are found important for the large mixing 
angle. We show that the new experimental data from BaBar provide a strong constraint on the 
mixing angle at the mass of Z)si(2536), from which two values can be extracted, i.e. 6i = 32.1° or 
92 = 38.4°. Our study agrees well with the latter one. Detailed analysis of the mass shift procedure 
due to the coupled channel effects is also presented. 

PACS numbers: 13.25.Ft, 14.40. Lb 



I. INTRODUCTION 



In the past few years one of the most important experimental progresses in the study of the charmed meson 
spectrum is the establishment of the lowest P-wave charmed-strange mesons, i.e. Z?so(2317), £'si(2460), I?si(2536), 
and Z3s2(2573) as now listed in Particle Data Group (PDG) 2010 Edition [T]. Since the first observation by BaBar 
Collaboration ,2j, the spin-0 state 1530(2317) and spin-1 Z3si(2460) (later confirmed by Belle [3] and CLEO ^) have 
initiated tremendous interests in its property and internal structure. These two states have masses lower than the 
potential model predictions, and their widths are rather narrow. It is somehow agreed that their low masses are 
caused by the open DK and D* K thresholds, respectively, and as a consequence, their narrow decay widths are due 
to the dominant isospin-violating decays, i.e. Z?so(2317) — >■ DgTT and 1^81(2460) — DItt (see the review of Refs. [11 16] 
and references therein). 

The heavy-light Qq system is an ideal platform for testing the internal constituent quark degrees of freedom. In 
the heavy quark limit the heavy quark spin is conserved and decoupled from the light quark degrees of freedom, 
which are characterized by the total angular momentum = -I- L, where is the light quark spin and L is its 
orbital angular momentum. With jq = 1/2 and jq = 3/2, one can arrange those four P-wave states into two classes, 
i.e. — 0+, 1+ and — 1+, 2+, respectively, where J is the meson spin as a sum of the heavy quark spin Sq 
and jq. For the axial vector states in the charmed and charmed-strange meson spectrum, since they are not charge 
conjugation eigenstates, state mixings between the ^Pi and ^Pi configurations are allowed. In the case of charmed 
and charmed-strange heavy-light system when the heavy quark symmetry is broken at order of l/rric, it would be 
interesting to study the mechanism that causes deviations from the ideal mixing scenario, i.e. breakdown of the heavy 
quark symmetry. This forms our motivation in this work. As mentioned earlier, Z3si(2460) and I?si(2536) lie near 
the threshold of D*K and both couple to D* K strongly via a relative S wave. It gives rise to coupled channel effects 
in the mass shifts of potential quark model calculations in comparison with the observed values and produces 

state mixings between the "^Pi and ^Pi configurations. Similar mechanism has been studied in the ao(980)-/o(980) 
mixing in Ref. jlO) . Determination of the mixing angle should be useful for understanding the property and internal 
structure of these two axial vector states. 
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We mention that various solutions have been proposed in the hterature to explain the observed results for Dsi(2460) 
and 1551(2536). For instance, D*K molecule or tetra-quark configuration have been investigated in Refs. pTHTS] , In 
Rcf. 14 , Z)si(2460) is explained as a dynamically generated state. The mixing angle has also been calculated in 
the quark model 1151 116] but with large uncertainties from the quark spin-orbital interactions. In this work, we 
investigate the two-state mixing propagator matrix which respects the unitarity constraint in a chiral quark model. 
We will show that the coupled channel effects via intermediate hadron loops can provide a simultaneous determination 
of the masses, widths and mixing angles of these two axial vector states. We also mention that the coupled channel 
effects on the ^Pi and ^Pi mixing was recently studied in Ref. |17) . where the the couplings were extracted in the 
^Pq model and a subtracted dispersion relation was applied to evaluate the hadron loops. In our approach we use the 
chiral quark model to extract the couplings and vertex form factor. We then extend the quark model form factor to 
a covariant form which can be applied on a general ground to much broader cases. 

The paper is organized as follows. In Sec. |lTj we give the basic formulas of two-state mixings via coupled channel 
propagators. In Sec. Ill the relevant coupling form factors are determined by the chiral quark model. In Sec. IV the 
propagator matrix is calculated in detail. Section |V] is devoted to show our numerical results for the mass and mixing 
parameters. The experimental constraints for the mixing angle are presented in Sec. |VI[ A summary is given in the 
last Section. In Appendix [X] the detailed definition and calculation of a special function used in the evaluation of the 
loop integrals with exponential form factors are provided. 



II. MIXING THROUGH COUPLED CHANNEL EFFECT 

We use \a) and \b) to present two pure states in the quark model. If they can couple to common final states, there 
will be a transition between them via single particle irreducible (IPI) diagrams as shown in Fig.[l] 
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FIG. 1: Transition through intermediate states 
The propagator matrix of \a) and \h) can be expressed as 



Ga,^[ ]S{\o)Ab)) 



The physical states \A) and \B) should be a mixture of \a) and 



IB) 



cose ~s\nee^^\( \a) \ _ ( \a) 

sin^e-^^ cose \\b) -^^'^'^> I \b) 



(1) 



(2) 



where R{e, (j)) is the mixing matrix, e is the mixing angle, and is a possible relative phase between \a) and \b). Then 
the propagator matrix of |^) and \B) is 



Gab — RGabR^ 



(3) 



The physical propagator matrix Gab should be a diagonal matrix. Thus, we can determine the mixing parameters 
{e, (j)} by diagonalizing the propagator matrix Gat- 
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In the present case, we set \a) = \^Pi), \b) = \^Pi), \A) = |I?,i(2460)) and \B) = |I?,i(2536)) as in Ref. The 
mixing scheme is 



|i:',i(2460)) = cose'l^Pi) -sin6'e^"^|iPi) 
|D,i(2536)) = sin6le-"^|3Pi) +cos6'|iPi) , 

where states pPi) and \^Pi) can be rotated to the eigenstates in the heavy quark hmit: 



2 /I 

3 V 3 

1 fl 

3 V 3 



The mixing angle 9 defined in Eq. Q can be related to 9' defined in j = 1/2 and j = 3/2 bases: 

e = 61' + 35.26° . 



(4) 



(5) 



(6) 



Considering parity conservation, the important intermediate states that can couple to i?si(2460) and Z)si(2536) 
are D*K, D*ri and DK*, of which the thresholds are fisted in Table |Tj 

TABLE I: The thresholds of intermediate states for i3si(2460) and Dsi(2536). 



Intermediate states 




D*+K° 




D"K*+ 




Threshold (GeV) 


2.501 


2.508 


2.660 


2.756 


2.761 



If all the particles involved are scalars or pseudoscalars. Fig. [T] will only represent sums of infinite geometric series 
and the resulting propagator matrix G becomes |10j 



Gab 



1 [ Db Dab 

DaDb - Dl, [ Dba Da 



(7) 



where Da and Di, are the denominators of the single propagators of \a) and respectively, and the mixing term 
Dab is the sum of all IPI diagrams, which satisfies Dba—Dab- But from Table |I| we find that the particles involved 
in the present case can be scalars, vectors or axial- vectors. There are five diagrams for the mixing of _Dsi(2460) and 
Dsi{2536) as shown in Fig. [2] 




FIG. 2: Mixing term for L»si(2460) and L>,i(2536). 
The mixing term can be generally divided into transverse and longitudinal terms: 



(8) 
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where P'"^ = g^'^ —p^p'^ /p^ and Q^^ = p^p^ /p^ are the transverse and longitudinal projector, respectively, and satisfy 

pM^^pA ^ pA^A^ Q^'^Qt = Q"^, P'^'Qt = . (9) 

Next we concentrate on the evaluation of the propagator matrix G^"^ for axial vector states. The numerator of the 
vector propagator is o^'' — p'^p'^/m^ and can be generally expressed as P**" + AQ^^ where A = 1 — p^ /m? . With 
the properties of Eq.([9]), the geometric sums, e.g. {a\S\h) in Fig. [l| can be taken for the transverse and longitudinal 
terms independently. After include the self-energy functions W^'^ and 11^'', the complete propagator matrix for the 
I'^Pi and l^Pi states becomes 

"^"^-^^ detG,,(.)+*'^ detG^, ' ^^^> 



with 



and 



Gabis) ^ Ml, - s^,s = < + - ^ , , (11) 

\ -Ilab{s) mi+na{s) - S I 



Gabis) ={ ' 2 I , (12) 



where M^^^ is the mass matrix. After diagonalization, the mass matrix becomes 

Ml^ = RMl,R^ ^ f ™^ \ \ . (13) 
Note that the longitudinal term G^^/detG^^ is nonvanishing, but the poles are only related to the transverse term 

Gab- 

By searching for the poles in the propagator matrix G^^{s), which is equivalent to set det G(s) — 0, we can obtain 
the masses and widths of the physical states. In general, there are two solutions sa and sb for the two state system. 
We can also extract the mixing angle 9a,b and the relative phase angle 4>a,b- These mixing parameters are different 
for these two states, since they are extracted at the physical masses of these two states, respectively. If G is a normal 
matrix, which means GG^ = G^G, then it can be diagonalized through a unitary transformation R. The resulting 
mixing angle 9 and relative phase (p can thus be uniquely determined. Otherwise, we can only get a quasi-diagonalized 
matrix through the unitary transformation R. The reason is because that orthogonality cannot be satisfied between 
these two physical states. 



III. COUPLING FORM FACTORS IN THE CHIRAL QUARK MODEL 

At hadronic level all the vertices in Fig. [2] involve the Axial- Vector-Pseudoscalar (AVP) type of coupling. In general, 
the AVP coupling vertex contains two coupling constants gs and go representing the S and D waves as shown in 
Fig.i 



\ 



FIG. 3: The AVP vertex via the S and D wave couplings. 
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Since the decay momentum is small near the threshold, we expect that contributions from the D-wecve coupling 
would be small. As a reasonable approximation, we omit go and keep gs to the order 0(u°). In the multipole 
approach, the helicity amplitude for 1+ — > 1^ + 0^ takes the form [T5] 



A, = {Sf,i^;0,0\S\S„iy) ^J2^L,0■,Sf,^I\S^,:y)YLaiq)GL 



(14) 



where Gl is the coupling constant for the L wave and q is the momentum direction of the final state particle in the 
center of mass frame of the initial state. In the present case, Eq. (14) becomes 



1 

/2n 
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In order to obtain gs to the order 0{v°), we set 

gs^Ao{q^O) = Aiiq^O) 



(15) 



(16) 



9j • q'j <ij • q'j 

c c c c 

(a) Emit a pseudoscalar (6) Emit a vector 



FIG. 4: Pseudoscalar (a) and vector meson (b) emission via an active light quark j in an effective chiral quark model. 



A. Coupling to D*K and D^r] 

One notices that at all the coupling vertices the interacting quarks involve only light quark, i.e. u, d and s. By 
treating the light mesons, pseudoscalar and vector mesons, as induced fields by a chiral Lagrangian for the mesons 
coupling to constituent quarks [TU], the light and heavy quark degrees of freedom can be separated out in terms 
of nonrelativistic expansions near the decay threshold. This approach has been successfully applied to light meson 
productions in photo-nucleon and meson-nucleon scatterings [20H26] and strong decays of heavy-light mesons [27j [28] 
recently. 

In the chiral quark model, we treat the pseudoscalar mesons K and 77 as the effective chiral fields as shown in 
Fig. |4][a). For emitting a pseudoscalar from an active quark line, the quark-meson coupling and corresponding non- 
relativistic form are respectively as follows [27^: 

ff,n = E ^4-V5.7^l7^V'ja^0™ , (17) 



3 



^™ ^ E 1^ {g^j- • q + hcT, ■ p; } exp(-zq • r,) , (18) 



with 
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where fm is the decay constant of the pseudoscalar meson, Ij the isospin operator, w the energy of the pseudoscalar, 
Mf and Ef the mass and energy of the final state heavy meson, a reduced mass given by 1/iJ.q = ^/nij + 1/m^-, 
Pj and Tj the internal momentum and coordinate for the light (jth) quark of the final state heavy meson. 

Following the procedure in [57], we derive the helicity amplitude A^, = {Sf,i'\Hm\Si,i') in the quark level. For 
I'^Pi — > l'^5'i + P, the explicit expressions are 



= igiha ex.p{- 



A? = z£ [2Gggi + M4a^ - g?)] exp(-^) 



and for l^Pi l^Si + P, we have 
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[2Gqqi + h{2a^ - qf)] exp(- 



(20) 



(21) 



where gi 



M/|/i |Mi) is the isospin factor, a the harmonic oscillator strength a = /3(2m2/(mi +1112))^^'^ as in 
Ref. [27], and qi = gm2/(mi + TO2). In the cs system, the 1st quark is s and the 2nd is c quark, and the flavor 
symmetry between the heavy and light quark is apparently broken. 

By taking equivalence between the quark and hadron level helicity amplitudes, we can extract the coupling form 
factor as follows: 



A 



. = ^{E, + M.){Ej + Mf}Al 



(22) 



Then from Eqs.(16), (20 1 and (21|, we finally obtain: 

6 

for l^Pi, 



gs = - — yZM,[l^f + Mf) ■ giAaexp(-^j , (23) 
for l^Pi, gs = j-^2Mi{Ef + Mj) ■ l=g,haeM~^) , (24) 
where i5 is a global parameter accounts for the strength of the quark-meson couplings as introduced in [27] . 



B. Coupling to DK* 



In this coupling, the vector meson K* is treated as an effective chiral field, for which the effective quark- vector-meson 
coupling Lagrangian and the corresponding non-relativistic coupling form [221 129] are 



H„ 



H. 



E 



2a*<? 



E 



i(Tj • q X e 



quj 



2m J Ef + Mf 2M'mj 



aij exp(— iq • r 



1/ ' 



qurri; 



/i ^ 2m J ' ^ 2Af> ^ 2M'nmj 



UJ 



2fiflq 



p) •q^a/jexp(-iq-rj) 



(25) 
(26) 
(27) 



where /i, uj and e are the mass, energy and polarization vector of the emitted vector meson, M' the sum of the 
constituent quark mass of the final meson, a the overall quark-vector-meson coupling, and other symbols have the 
same meaning as those in Eqs.( 17 1-( 19 ). Using the above operators, we can extract the helicity amplitudes AI, i.e. 
for l^Pi ^l^So+W, 



A\ 



^agiB— exp{--—r), A^ = 

2y/2 a ^ 4a2^ ° 



and for l^Pi ^ l^^o 



Al ^ iagiAacxp{-—^), A^ = —agi 



-C—+Da 



2a2 



- 1 



(28) 



(29) 
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with 



A 



B 



1 



c = 

D = 



^2q 

q 



1 

2 m,, 



Ef + Mf 



-) + 



quo 



2nHq 



2M'm 



qujm 



3 

I -I 



2M'fi 2M'fimj 



Substituting Eqs. ([281 and ^ to Eqs. ([16|) and ([221 , we obtain 
for l^Pi, 55 = 0, 

for fiPi, gs = -J{E,+M,){Ef+mf) 



1 a , 
agi — exp(- 

fJ-q 



V2 



4q2 



) 



(30) 

(31) 
(32) 
(33) 

(34) 
(35) 



The coupUng gs — ^ hi Eq. (34) is because A\ in Eqs. (28) is proportional to qi. Thus, the effective couphng vanishes 
below the open decay threshold. As a consequence, the contributions from Fig. [2](4) and (5) should vanish to the 
order 0(ti°). Therefore, we only need to consider the contributions from Fig.[2]jl-3) in the following calculations. 



C. Numerical results for couplings 



In the numerical calculation, we set \(f\ =0 when the initial state lies below the threshold for VP final state [29j . 
We adopt 77 = ^[uu + dd— ss) in the ij — ij' mixing scheme which corresponds to 9p = — arcsin(l/3) — —19.47° for 
the flavor octet and singlet mixing. Since the contribution from the D*?] loop is small, the uncertainties with 9p have 
only negligible effects on the mixing matrix. We obtain isospin factors gi for different intermediate states as listed in 
Table nil 



TABLE II: Isospin factors gi extracted in the quark model. 



I^PI/I^PI 


D*'>K+ 
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1 


1 


1 

73 



The following values are adopted for other parameters [57]: 6 — 0.557, (5 = 0.4 GeV, Jk = f-q = 160 MeV, and 
the constituent quark masses m„ = md — 350 MeV, = 550 MeV, nic — 1700 MeV. We note that our numerical 
results are not sensitive to nic = 1500 ^ 1700 MeV, while the light quark masses m„ — md = 330 ^ 350 MeV and 
ms = 500 ~ 550 MeV will lead to about (1 ^ 5)% uncertainties with the final results. 

The masses Mi of the initial states l^Pi and l^Pi cs still have uncertainties. Fortunately, the couplings \gs\ change 
only 5% at most when Mi € [2.460, 2.536] GeV as shown in Fig. [s] Also it shows that the couplings to D*^K~^ and 
are almost the same for each state due to the isospin symmetry. A set of typical gs couplings is listed in 
Table mi 



TABLE III: Vertex couplings gs at Mi = 2.5 GeV. 



gs (GeV) 








l^Pi 


-7.982 


-8.052 


2.040 


l^Pi 


5.644 


5.694 


-1.443 



Apart from the on-shell coupling 175, the chiral qua rk mo del als o p rovides an exponential momentum-dependent 
form factor exp{—q1/Aa'^) as shown in Eqs. (23), (24), (34|, and (35). In order to keep this feature in the meson 



loops, we modify the exponential form factor to a covariant form: 



9 99 
I it \ . ,q~m\ 



(36) 
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FIG. 5: the absolute values of couplings gs as functions of the initial meson mass Mi 



where q'^ and m are the four-vector momentum and mass of either V or P particle. Parameter A is the cut-off energy, 
which can be determined by the quark model, namely, for the D* K and D*rj loops, 

/? = 1.174 GeV. (37) 
The exponential form factor serves to remove the ultraviolet divergence in the loop integrals. 



A = 2- 



2m„ 



IV. THE PROPAGATOR MATRIX 



In this Section we will determine the propagator matrix G. From Eqs. (23) and (24), we have 

n, = -Vmab, lib = -nab/V2. 



(38) 

So, we only need to calculate the mixing term Ilab. With the AVP coupling form factors, we can explicitly write down 
Z?^^ as the following: 



D^b 

ab 



X! 9agb ■ 



d^k exp(^^^j exp(^ 
(27r)4 [k^ - ml][{k + pY ^ ml 



A2 



= J2 gagbiHg^" + B'p'^p'^/p^) = ^ gagbloop , (39) 
where ga and gb are the iS-wave couplings of the two vertices, respectively. Comparing with Eq. (|8]), we obtain 

llab = J2gagbU . (40) 
The mixing term Ilab can be decomposed into two terms, i.e. 



n 



ab) 



with 



(41) 
(42) 



where 11 and 11^^ are contributions from the g'"^ and fc^fc" terms of the vector propagator, respectively. 

As follows, we first make an on-shcU approximation to investigate the absorptive part. Then, we investigate the 
full integrals with the help of the exponential form factors. 
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A. On-shell approximation 



Since the absorptive part of a two-point function is independent of the form factors, the on-shell approximation 
will allow us to separate out the absorptive par t an d then compare it with that in a full loop integral. Here we only 
consider 11^ , for which the loop integral of Eq. ( 39 1 in the on-shell approximation becomes 



loop, g^^-l-l^BO{s^ ml, ml) ^ U^g^"^ 



The resulting mixing term Ilab is a function of s. We plot Tlabi\/s) in Fig- |6]with the couplings listed in Table 
adopted. 



(43) 

m 



> 

O 




2.60 



FIG. 6: The mixing term Hat in the on-shell approximation. 

In Fig. |6j two kink structures can be identified. The first one at y/s — 2.501 GeV corresponds to the D*^K^ 
threshold, and the second one at y/s — 2.508 GeV to the D*^K'^ threshold. This result will be compared with the 
absorptive part in the full loop integrals later. 



B. Full loop calculation with the exponential form factor 



In this Subsection we perform the full loop calculation with the exponential form factor. The explicit formula for 

ni is 



loop I 



d^k exp(^^r^j exp(^ 
(27r)4 [fc2 -m2][(A:+p)2 -to2] 



A2 



with 



c 
A 



d2:e"L/(2,l,— ,aA) =11^^ 
a 



2 

A2 ' 

(1 - 2x)2 

s(2.T^ — 2a; -I- 1) — — m^ 
(1 — x)m1 + xmp — x{l — x)s 



(44) 
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The explicit formula for 11^ is 



loop2 = (-i) 



(^) 



exp 



( jk+p) 



A2 



(2^)4 [fc2„TO2]p-+p)2-TO2 



(45) 



where a, 6, c, A are the same as those in Eq. (44). The function [/(a, 6, c, z) is a class of special integrals which 
appears in the evaluation of the loop integrals with exponential form factors. The detailed definition and calculation 
of U{a, b, c, z) are provided in Appendix [A| 

The full loop calculation of the mixing term Tlab{\/s) is presented in Fig. [7] where the parameters are the same as 
before. In order to see clearly the contributions from different parts, we also give two sets of the calculated values in 
Tables [IV] and IVl 
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FIG. 7: (color online). The mixing term with exponential form factors. 11^;, and 11^^ are the contributions from the 
g^^ term and k^k" term of the vector propagator, respectively. The dashed lines represent the dispersive parts, while the 
dot-dashed lines represent the absorptive ones. 



TABLE IV: The mixing term Eat at the pole position = 2.4545 GeV. 



intermediate state 






Din 


na5 = E<7.3i>n (GeV^) 


5.5i,(GeV2) 


-45.05 


-53.90 


-2.944 




n^(on-shell) 
















-2.391 X 10"^ 


-2.253 X 10"^ 


-0.868 X 10-=* 


0.2317 




5.650 X 10"^ 


5.448 X 10"^ 


2.503 X IQ-'^ 


-0.0056 



TABLE V: The mixing term n^i, at the pole position = (2.5449 - O.OOlOz) GeV. 



intermediate state 








nab^E^afft-n (GeV') 


.gaff6(GeV2) 


-45.05 


-53.90 


-2.944 




(on-shell) 


-2.969i X 10"^ 


-2.718i X 10~^ 





0.2803i 




(-4.157- 2.969i) x 10~^ 


(-4.260- 2. 718i) x 10"^ 


-1.389 X 10"^ 


0.4210 + 0.2803i 




(9.419 + 0.886i) x 10"^ 


(9.059 + 0.678i) x lO"'^ 


3.560 X 10"^ 


-0.0092 - 0.0008i 



The loop calculation results help us to learn the following points: 
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• The imaginary part of 11;^^ with exponential form factors is the same as that in the on-sheh approximation. It 
justifies our calculation method for U{a, h, c, z) as described in the Appendix. 

• The contribution from the term oi g^'^ is dominant. The open thresholds of D*^K'^ and D*^K^ cause two kinks 
in both real and imaginary parts. With the increase of ^/s, Ren;^;^ first increases until it reaches a summit at the 

threshold. It then decreases in a linear behavior in terms of -^s. In contrast, Imll^^ is zero below the 
£)*o^+ threshold and then increases quickly when the decay thresholds are open. One can see that below the 
D*^K^ threshold, the real part is the only contribution and cannot be neglected. The imaginary part becomes 
significant above 2.53 GeV. 

• The calculation also shows that the contributions from the k^k'^ term of the vector propagator are negligible. 
Near the threshold, the momentum is small such that 11^^ suffers an 0(l/m^) suppression comparing to n;'^;^ in 
both the absorptive and dispersive part. 

• The contributions from the D* K loops are found dominant, while the contributions from D*i^ account for only 
about 1% of the mixing term due to the rather small coupling value in the D*7] loop. 



V. POLE POSITIONS AND MIXING PARAMETERS 



With the Iiab{s) determined, we can directly search for poles for the physical states in the propagator matrix G in 
Eq. (10). We adopt the following bare cs masses, m^Pi] = 2.57 GeV and m[^Pi] = 2.53 GeV, from the Godfrey-Isgur 
(GI) model HHj as input. By scanning over the energy y's, the requirement of | det[G(s)]| = provides a direct access 
to the pole positions as shown in Fig. [8} T wo possible poles near 2.46 GeV and 2.54 GeV are highlighted. When 
varying the cut-off parameter A in Eq. (371 within the range of [1.174 — 0.22, 1.174 + 0.22] GeV , it shows that the 



higher pole is stable and the lower one changes from 2.47 GeV to 2.44 GeV. Searching for the poles on the complex 
energy plane, we can pin down the masses and widths of these two poles as listed in Table [Vl] It shows that the mass 
of Z?si(2460) changes 3.6 MeV at most with or without the contribution from the fc^fc'' term of the propagator, while 
the mass of Z?si(2536) changes only 0.1 MeV. The extracted mass of Dsi{2AQQ) is only 5 MeV below the experiment 
value, and the mass of Dsi(2536) is only 10 MeV above the experiment one. In principle, the Okubo-Zweig-Iizuka 
(OZI) rule allowed hadronic decay width of Dsi(2460) is zero. The obtained width 2.0 MeV for Dgi{2b'iQ) seems to 
be slightly larger than the experiment value 0.92 MeV, but can still be regarded as in good agreement. In brief, our 
prediction for the masses and widths of these two states agrees well with the experiment data. 




2.40 



2.45 



2.50 2.55 

rr(GeV) 



2.60 



FIG. 8: Pole structures highlighted by the zero values of det[G'] in the propagator matrix. 



Before extracting the mixing parameters, we show that our formalisms can reproduce the ideal mixing angle in 
the heavy quark limit. In this limit, and m;, are degenerate. From Eqs. (11 1 and (381, we only need to diagonalize 
the simple matrix 




(46) 
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TABLE VI: Masses and widths obtained from the pole analysis. 



[m-i^] (MeV) 


Dsi(2460) 


_D,i(2536) 




2454.5 


2544.9 - l.Oi 




2455.8 


2544.9 - l.li 


Experiment 


2459.5 


2535.08 - 0.46i 



which immediately leads to 9q = arctan[l/\/2] = 35.26°. 

Now we proceed to the extraction of the mixing parameters {6,4>} by diagonalizing G{s) with ^/s fixed at the 
poles. When G is a complex matrix, we try to approach the diagonal limit RGabR^ = Gab in three ways: Method 
I, set Gfa^ = 0; Method II, set G^f = 0; and Method III, minimize IGfj^l + \G^f\. The results from these 
three diagonalization schemes are listed in Table |VII| As we expected before, the mixing angles of these two states 
determined at their pole masses are indeed different. For Dsi(2460), G is a symmetric real matrix. So the mixing 
parameters are the same in these three methods: 9 ~ 47.6°, </> = 0°. From the mixing scheme in Eq. Q, 61 > 45° 
means that the ^Pi component is larger than the ^Pi in _D5i(2460). This mixing pattern would affect the mass shift 
as we will show later. The result corresponds to 9' — 12.3° in the j ^ 1/2 and j = 3/2 mixing in the heavy quark 
limit. For Dsi{253Q) , G is a complex matrix. The mixing angle 9 = 39.7° determined at the Z3si(2536) mass changes 
little in those three methods, while the relative phase suffers an uncertainty of = —6.9° ^ 6.9°. We will show 
later in Sec. |VI| that the mixing angle 9 — 39.7° is consistent with the experimental constraints and can be useful for 
picking up one of those two solutions from the experimental fit. Again from the mixing scheme, 9 < 45° means that 
the ^Pi component is larger than the '^Pi one in Dsi(2536). The result corresponds to 9' = 4.4° in the j = 1/2 and 
j — 3/2 mixing bases. The energy dependence of the mixing angle reflects the breaking of orthogonality among these 
two physical states. 



TABLE VII: The mixing angle 9 and relative phase extracted at the two poles in those three diagonalization schemes. 





Dsl(2460) 


I>si(2536) 


I 


II 


III 


I 


II 


III 




{47.5, 0} 


{47.5, 0} 


{47.5, 0} 


{39.7, -6.4} 


{39.7, 6.4} 


{39.7, 0} 


nib + Ilab 


{47.6, 0} 


{47.6, 0} 


{47.6, 0} 


{39.8, -6.5} 


{39.8, 6.5} 


{39.7, 0} 



From the mixing angle analysis, we also learn that the Dsi{24:60) has a larger j — 1/2 component which couples to 
the D*K through an S'-wave. It hence acquires a significant mass shift ^ 100 MeV through meson loop corrections. 
In contrast, the I?si(2536) contains a larger j = 3/2 component which couples to the D*K through a ZJ-wave. It only 
gains a small mass shift ~ 10 MeV. 
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Z),i(2460) 
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FIG. 9: (color online). Schematic plot for the mass-shift procedure. The thin solid bars represent the original ^Pi and '^Pi 
states in the quark model. The thick solid bars represent the two physical states Dsi(2460) (left) and Dsi{25S6) (right). The 
solid arrows represent the mass shifts due to the diagonal elements Ha and Hb, while the dashed arrows represent those due to 
the off-diagonal element Hat- The threshold for D* K is shown by the horizontal dashed line. 



The mass shift procedure is also an interesti ng issu e and can help us to understand why P)si(2460) has a larger 
Pi component. As shown in Fig. |9] and Table VIII we can decompose the mass shift procedure into two classes. 
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TABLE VIII: Mass shift procedure at different ^/s. From we can see the mass shifts due to Tla and Tib, while from 

\g , further mass shifts due to Ilai, can be learned. Here we only present the results from the quasi-diagonalization method 
II. 



^ GeV bare mass (GeV) -/M^ (GeV) \/A^1b (GeV) 



2.4545 GeV 








2.497 


0.481i 






2.548 





















2.53 







0.481i 


2.505 









2.455 




2.5449 GeV 





2.57 




2.471 - 


- 0.040i 


0.206 - 


0.678i 




2.545 - 


O.OOli 


0.141 - 


h 0.303i 










0.206 - 


- 0.678i 


2.454 - 


O.OSOi 







2.379 - 


- 0.123i 





i.e. diagonal shift and off-diagonal shift. The diagonal elements Ha and lib cause both ^Pi and -^Pi states to 
move downwards, while the off-diagonal elements Hab make one state to shift up and the other to shift down. At 
= 2.46 GeV, after the diagonal shift the '^Pi state is still higher than the ^Pi. But after the off-diagonal shift, 
the higher mass state moves down to become an on-shell £'si(2460) and the lower state moves up to become a virtual 
Dsi(2536). The reversal of the mass ordering results in a mixing angle 9 > 45° and thus a larger ^Pi component 
in Dsi{2A60). At y/s = 2.54 GeV, after the diagonal shift the ^Pi becomes higher than the "^Pi. Then after the 
off-diagonal shift, the higher state becomes much higher and the lower much lower, which causes a mixing angle 
6 < 45° and a larger ^Pi component in Dsi(2536). Note that in this situation the on-shell state corresponds to the 
1551(2536), and the Dsi(2460) appear as a virtual one. 



VI. EXPERIMENTAL CONSTRAINTS ON THE MIXING ANGLE 



In this part, we come to survey the constraints for the mixing angle 6 from experiments. The strong decays of 
Dsi(2536) has been measured with reasonable precision which are summarized in Table IX Since the D*K channel 



is the only allowed strong decay channel for £'si(2536), it is a good approximation to assume 

r[i:>,i(2536)] « r(i:»3i(2536) ^ d*k) , 



(47) 



which can be estimated in the chiral quark model. The partial width fractions Pi and R2 can also be calculated and 
compared with the data. 



TABLE IX: The available experimental status of _Dsi(2536). 



m = 2535.0 


8 ±0.01 ±0.15 MeV, 


r = 0.92 ±0.03 ±0.04 MeV (BaBar g^) 




n _ r(I5*(2007)"if+) 
^ r(D*(2010) + X0) 


= 1.36 ±0.20 (PDG2010 [T]) 


R2 


r(D* (2010)+ K")s-^ai 
^ r(D*(2010) + KO) 


' - 0.72 ± 0.05 ± 0.01 (Belle EB) 



The hehcity amplitudes for l^Pi — > D*K and l^Pi D* K have been hsted in Eqs. (20) and (21). The partial 
width can be obtained by [27] 



r = 



AttM, 



(2J,±1)^' ' 



(48) 



where Ji is the spin of the initial particle. In order to calculate P2, we need to extract the 5- wave components from 
the helicity amplitudes. By defining As — and Ad — we deduce from Eq. (151 



^0 = - V2Ad 



^1 = ± -\Ad 



As = \{A^ + 2Ai) 
Ad = -^{Ao-Ai) 



(49) 



where the S and D-wave components have been separated out. We use the same model parameters as before to 
calculate the partial width r[I?si(2536)] and ratios Pi and P2 in terms of the mixing angle 9. The results are shown 
in Figs. ([10])-([13]). 

A similar result as Fig. 10 for r[£'si(2536)] in terms of 9 has been given in Ref. [27] but with the notation 
6' — + 90°. Those three horizontal lines in Figs. (10)-(12) represent the upper limits, center values, and lower 
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FIG. 10: r[D,i(2536)] as a function of 9 FIG. 11; Ri as a function of i 



FIG. 12: i?2 as a function of 6* 



limits of the experimental data. The interesting feature arising from the results of Figs. ( 10 )-( 12 1 is that the overlaps 
between the experimental data and theoretical values are separated into two narrow bands of 9 which are located 
symmetric to the ideal mixing angle Oq = 35.26°, i.e. 9i ~ 32.1° or 6*2 — 38.4°. An alternative way to present the 
results is via Fig. |13[ where the overlapped 6 values are denoted by the vertical dashed lines, while the experimental 
observables with errors are presented in terms of 9. Notice that these two bands of 9 are both smaller than 45°. 
Therefore, based on the present experimental measurements, one cannot determine which value for 9 should be taken. 
It turns out that our analysis in Sec. |V] can precisely pick up one of these two solutions, namely, 02 — 38.4° is favored 
in comparison with the theoretical value 9 = 39.7°. 



Ri 



28 30 32 34 36 38 
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FIG. 13: The experimental constraints on the mixing angle 9. 



VII. SUMMARY 



In summary, we have studied the mixing mechanism for the axial vector states 1581(2460) and 1581(2536) via the 
S'-wave intermediate meson loops. We establish the propagator matrix for this two-state system. Then, by searching 
for the pole structures in the propagator matrix, we can pin down the masses and widths of the physical states. 
The mixing angle and relative phase between the '^Pi and ^Pi components can be determined by diagonalizing the 
propagator matrix. For i:»si(2460), we obtain m = 2454.5 MeV, 9 = 47.5° and (j) = 0°. For Dsi{2536), we find 
TO = 2544.9— l.Oi MeV, 9 = 39.7°, and (j> = —6.9° ~ 6.9°. Our results agree well with the experimental measurement. 
In particular, the new BaBar measurement put a strong constraint on the mixing angle at the mass of I?si(2536) with 
two solutions, 9i ~ 32.1° and 02 — 38.4°. Our theoretical calculation finds 9 = 39.7° which is in good agreement with 

02- 

Note that due to the breaking of orthogonality the energy-dependent mixing angles defined at the different physical 
masses turn out to have different values. We find that both iDsi(2460) and Z)si(2536) have a relatively large ^Pi 
component in their wavefunctions. 

It is also interesting to learn the important role played by the coupled channel effects for states near open thresholds. 
For states that can couple to each other via the coupled channels, the two-state propagator matrix carries rich 
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information about the mixing and mass shifts as a manifestation of the underlying dynamics. Extension of such a 
study to other axial- vector meson mixings would be useful for deepen our understanding of the coupled channel effects 
and their impact on the hadron spectrum. 
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Appendix A: Calculation of function U{a,b, c, z) 

Initially we define 

U{a,b,c,z)^^^l^^ dtt--\l + tf-'^''cxp(^-zt-^^ , (Al) 

which is the typical integral we encounter in the calculation. A special case, U{a,b,c — 0,z) — U{a,b,z), is the 
Tricomi confluent hypergeometric function, which is a build-in function in Mathematica. The function U{a, b, z) is a 
single- valued function on the z-plane cut along the interval (— cxi,0], where it is continuous from above, i.e. 

when z < 0, C/(a, b, z) = lim Uia, b,z + ie) . (A2) 

Function U{a, b, c, z) as a physical quantity should be analytic with respect to its arguments. However, the integral in 
Eq. (Al) only converges when Re{z) > and Re{a) > 0. In o rder to analytically continue the integral to Re{z) < 0, 
we make a change in variables zt = x{z > 0). Hence, Eq. (Al) becomes 

1 — poo / 
Z I , „_^ , xfc-Q-1 / CZ 



[/(a, 6, c, z) = ^ y da; x''-\x + z)""""^ exp ^-x - ) . (A3) 



In the region Re{z) > 0, Eq. ( |A1[ ) and (A3) are exactly equivalent to each other when Re{b) < 2. The difference 
between them can be expressed by the integral over Cr in Fig. [T4|(a). 




FIG. 14: The continuation of U(a,b,c,z) 



(All 



the integral in Eq. ( A3 ) has larger 

-z 



When Re{h) < 2, the contribution from Cf> is zero. Comparing with Eq 

convergent region, i.e. the whole complex z-plane except z < 0. When z < 0, there is a singular point at x 
in the integral path as shown in Fig. |14[ b). Considering Eq. (A2), U{a,b,c, z) must satisfy a similar requirement. 
It means that the integral path Ci in Fi g. [T4| (b) should be replaced by the integral path C2. Using the expression 
in Eq. (A3) and the replacement in Fig. |14[b), we can analytically continue the integral in Eq. (All to the whole 
z-plane. The constraints of the above method are Re{a) > and Re{b) < 2, which could satisfy our present need. 
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To test this method, we compare the results for U (a, b, z) in Fig. 15 using our method and the build-in Mathematica 



program. It shows that these two calculations are in good agreement to each other. This test is done at c = 0. Since 
c in Eq. (A3) does not bring either new divergence problems or new singular points, we can justify that our analytic 



continuation of C/(a, 5, c, z) is quite reliable and generally applicable. 
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FIG. 15: Test the accuracy and precision of the analytic continuation of U{a, b, c, z) 
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